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Abstract 

Our aim is to compute the lower moments of the unpolarized and polarized deep- 
inelastic structure functions of the nucleon on the lattice. The theoretical basis of 
the calculation is the operator product expansion. To construct operators with the 
appropriate continuum behavior out of the bare lattice operators one must absorb the 
effects of momentum scales far greater than any physical scale into a renormalization of 
the operators. In this work we compute the renormalization constants of all bilinear 
quark operators of leading twist and spin up to four. The calculation is done for 
Wilson fermions and in the quenched approximation where dynamical quark loops are 
neglected. 
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1 Introduction 



This work is part of an ongoing effort [f , 2] to compute the deep-inelastic structure func- 
tions of the nucleon both for unpolarized and polarized beams and targets. The theoretical 
basis of the calculation is the operator product expansion (OPE). The OPE relates the 
moments of the structure functions to forward nucleon matrix elements of local operators. 
Lattice simulations of these matrix elements, combined with a perturbative evaluation of 
the short-distance parameters of the expansion, i.e. the Wilson coefficients, can provide 
complete information of the deep-inelastic structure of the nucleon. 

A necessary ingredient of the calculation is the renormalization of the lattice operators. 
The bare lattice operators are in general ultra-violet divergent. They must be renormalized 
such that in the limit of zero lattice spacing the renormalized operators correspond to finite, 
Lorentz covariant operators which obey the same renormalization conditions as those in 
the continuum. In particular, operator matrix elements and Wilson coefficients must be 
computed in the same renormalization scheme. 

In principle the renormalization constants can be computed in perturbation theory, like 
the Wilson coefficients. However, the lattice operators may be subjected to large tadpole- 
induced renormalizations which lead to poor convergence of the perturbative series, and 
the bare coupling constant may not be the most effective expansion parameter to use [3]. 
Both problems can be remedied by simple redefinitions of the basic operators used to define 
the lattice theory and the expansion parameter, once the perturbative result is known. 
The validity of this procedure can finally be checked by computing the renormalization 
constants non-perturbatively [4, 5]. 

In this paper we study the renormalization of quark bilinear operators in perturbation 
theory to one loop order. The renormalization of gluonic operators, which enter the OPE 
to order a 2 s , will be dealt with in a future publication. We shall focus on operators of 
leading twist and spin less or equal to four. We use Wilson fermions and work in the 
quenched approximation where one neglects the effects of dynamical quark loops. A brief 
account of our results was given in ref. [2]. 

The paper is organized as follows. In section 2 we present the operators and discuss 
their transformation properties under the hypercubic group. In section 3 we state the 
renormalization conditions. Section 4 contains a brief description of the calculation. The 
results are presented in section 5 for a few representations which are particularly suited 
for numerical simulations. Finally, in the appendix we give results for the full tensors, 
from which the renormalization constants for all other representations can be deduced. 
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2 Operators and Representations 



The cross section for deep-inelastic lepton-nucleon scattering can be written in terms of 
four structure functions, F\, F 2 , g\ and g 2 . Here F\,F 2 are the well-known spin averaged 
structure functions and gi,g 2 are the polarized structure functions. 
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where the superscript / distinguishes the various operators and its Wilson coefficients, 
and fj, denotes the subtraction point. In eqs. (2) - (4) {• • •} ([• • •]) means symmetrization 

(anti-symmetrization) of the indices. For the definition of C^^]..^ } see re f- [6]. In the 
following we shall restrict ourselves to the leading quark bilinear operators 



Oln-nn = 9 ^1<j 75 f w ••• D„ n if} - traces. (6) 



We may do so because we will be working in the quenched approximation. The corre- 
sponding gluon operators are 

Ou,...u n = l n ~ 2 TrF"D, 2 ---D, n _,F a , -traces, (7) 



— ' - L1J -fl 1 - Ly ^2 ^ fln-l-L a/In 

Ot, „ = i^TrF" D u , ■ ■ ■ D u ,F au -traces, (8) 
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where = \t llvpa F p ' J . We will deal with the gluon operators in a future publication. 

The traceless and symmetric operators C^^...^} and O^^...^ y transform irreducibly 
under the Lorentz group. The r.h.s. of eqs. (2) and (3) are the only traceless, symmetric 
tensors of maximum spin, n and n+ 1, respectively, one can build from a single momentum 
vector and the polarization vector s^. Both operators have twist two. The operator 
^[alm]---!! }' wri ich is also traceless but of mixed symmetry, transforms irreducibly as well. 
Again it is only possible to form a single traceless, mixed symmetric tensor from the 
momentum and polarization vectors, as given by the r.h.s. of eq. (4). This operator has 
spin n and twist three. 

The moments (1) are independent of fj,. The operator matrix elements inherit the 
//-dependence from the renormalization constants, which we will introduce in the next 
section. This means that the //-dependence of the renormalization constants must 'match' 
the //-dependence of the Wilson coefficients c^|(// 2 /Q 2 , g(fi)), e^|(// 2 /Q 2 , g(fi)) which are 
known perturbatively. 



Lattice Operators 



We do a Wick rotation to obtain the operators in euclidean space-time. On the lattice we 
furthermore have to replace the covariant derivatives by the lattice covariant derivative 

D^(x) = ^ (U x ^(x + A) " Ul_^(x - //)) (9) 
^(x)D^ = (Hx + fiWlf,- i){x - {i)U x -p,^ , (10) 
where a is the lattice spacing. The link matrix U Xtfl is given by 

U Xltl = expiiagA^x)}, = T a A a ^ (11) 
where g is the bare gauge coupling and T a are the generators of the SU(3) Lie algebra. 
We use the Fourier decomposition of the fermion fields and the gauge potentials 




where the momenta are restricted to the first Brillouin zone. The phase factor e M ?W 2 i n 
eq. (14) amounts to redefining the potentials at the midpoints of the links connecting two 
neighboring lattice sites. In momentum space the operators then take the following form: 



4 



Order g° 

n ■ / \ 

o^ n = ^(p)t w ^(p) n sm - traces ; ( 15 ) 

i = 2 a 

Order g 1 



i = 2 3 = 2 k=i+l 

X A w (?)cos(a(p' + p w )/2) - traces, (16) 



where p' = p + q; 
Order g 2 



( n-l n i-1 j-1 n ^ m ( nn ' s 

o^ n = g*m^m\ EEnn n ( ^ 



k i=2 j=i+l k=2 l=i+l m=j + l tt 

v sin(a(p w +g w ))sin(ap w ) 

a a 

X cos(a(p' w - g ;./2)) ^.(g) cos(a( PAlj + (17) 

2 8 = 2 j = 2 k=i+l a a 

X A^) sin(a(p' w +P«)/ 2 )| " traces, 

where p' = p + q + q'. 

The operators O^^...^ are obtained from eqs. (15) - (17) by replacing V'(p')7^iV'(p) 
by ipijj'Yj^j^ipijj) and by adjusting the rest of the indices appropriately. 

These expressions for the operators, together with the well-known lattice Feynman 
rules [7, 8], will form the basis of our calculation. 



Representations 

In euclidean space-time the Lorentz group is replaced by the orthogonal group 0(4), 
which on the lattice reduces to the hypercubic group H(A) C 0(4). Accordingly, the 
lattice operators are classified by their transformation properties under the hypercubic 
group and, of course, charge conjugation. 

In the continuum it was relatively easy to identify the operators of leading spin trans- 
forming irreducibly under the Lorentz group. The corresponding lattice operators will 
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in general not transform irreducibly, which allows them to mix with lower- dimensional 
operators under renormalization. A necessary condition for a lattice operator to be mul- 
tiplicatively renormalizable is that it belongs to an irreducible representation of H(4). 

Under charge conjugation the operators transform as 

C/i!/i 2 ---/in ~> ( _ 1) C/il/in/in-l-"/"2 , (18) 

o 5 — ► (—^) n o 5 (^ ( )) 

Usually only C = + operators contribute to the OPE (1). In the quenched approximation, 
however, also operators O l _ tl ... l _ tn with odd n are relevant. 

In ref. [9] we have identified all irreducible representations carried by the operators 
O and O 5 up to rank four. Furthermore we have given explicit bases for all irreducible 
subspaces. We will summarize our results here. 



We follow the notation of ref. [10] (see also ref. [11]) and denote the irreducible represen- 
tations of H(4) by rj!\ where / is the dimension of the representation and k distinguishes 
inequivalent representations of the same dimension. The defining representation will be 
called t[ 4 \ Let us first consider the operators O^...^, which transform as the tensor 
product 



(gK (4) - 
i 

For the traceless tensors of n = 2, 3 and 4 we obtain the decomposition 



(20) 



r< 4 > ® r< 4 > C = + rf>8rj 6 >8r< 6 >, (21) 

S v ' 

symmetric 

r< 4 > ® r< 4 > ® r< 4 > C = r< 4 > 8 r 2 (4) 8 r< 8 > 8r< 4 > 8 r< 8 > 8 rf\ (22) 

S v ' 

symmetric 

r| 4) 8 r< 4 > 8 r< 4 > 8 r< 4 > C = + r' 1 * 8 r' 1 * 8 r< 2 > 8 r< 3 > 8 r< 6 > 8 r< 6 > 8 r< 6 > 8 r< 2 > 8 r< 2 > 

S v ' 

symmetric (23) 
82r 2 (3) 8 27-f 8 rf 8 2r 1 (6) 8 3r 2 (6) 8 2^ 8 ^ . 

We have distinguished between representations on totally symmetric tensors and tensors 
of different symmetry. We see that, unlike in the continuum, the representations are not 
uniquely determined by requiring, e.g., that the tensors are traceless and symmetric. There 
are several representations which appear twice, e.g., in eq. (22), and we will see that 
they can mix due to the presence of non-0(4) covariant contributions, which complicates 
the matter considerably. 

Let us now consider the operators Ol lll ... /ln and their representations. Here one factor 
of gets replaced by ® rj 1 ^ = rj 4 \ Thus the operators transform as the tensor 
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product 

n+l 

{®t1 4) )®tP. (24) 
1 

We are interested in the cases n = and 2. For the traceless tensors we find the decom- 
position 

rl^rf) (25) 
r< 4 > ® r< 4 > ® r< 4 > ® r' 1 ) rj 4 » ® r< 4 > ® r< 8 > 8r< 8 > ® r< 8 >. (26) 

S v ' 

symmetric 

Again, there is the possibility of mixing under renormalization. Unlike in the continuum, 
operators of the form (5) may also mix with operators of the form (6). 

The numerical simulations are greatly facilitated if one keeps the spatial momenta as 
small as possible. We shall therefore concentrate on those representations which require 
a non-zero external momentum in at most one spatial direction. The renormalization 
constants of the operators in the other representations will be given in the appendix. 

We assume that the spatial momentum is either zero, or non-zero only in the 1- 
direction. Following ref. [9], we are then led to consider the operators and representations 
listed in table 1. It will turn out that it is not sufficient to consider the totally symmetric 
operator 

C{114} - 2(^{224} + C{334}) (27) 

because it mixes with the mixed symmetric operator 

«4n» ~ 2 < ^°((' i22 )) + C '{{ 433 »)' ( 28 ) 
as indicated by the curly brackets, where 

@{{fivp)) ~ @p,vp "I" & p,pv ~ Opp,v ~ Qpvp,- (29) 
(8) 

Both operators belong to the representation r| In case of the operator O^^y we will 
consider two different representations. In the continuum limit both of them must lead 
to the same final result. The purpose is to check for finite cut-off effects in combination 
with our numerical simulations. In the last column of the table we have listed the matrix 
elements to which the operators contribute. Note that we have changed our notation. The 
operator (28) is the negative of the operator considered previously [2]. 



3 Renormalization Conditions 



Let us denote the lattice regularized operators by 0(a). We then define finite operators 
0(fj,), renormalized at the scale /i, by 

O l (ji) = Z l ^((a^)\g(a))O m (a), (30) 
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Rank 


Operator 


Representation 
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(O) 


2 
2 

3 

4 

1 

3 
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C{14} 

C{44} - K {H} + ^{22} + C>{33}) 

| C{114} - \(0{224} + C{334}) ] 
1 °«411» " |(°«422» + C«433») J 

C{1144} + C{2233} - C{1133} ~ C{2244} 

U 2 
(I 5 

U {214} 
U [2{1]4} 


r (6) 
r 3 

(3) 

(8) 
T l 

(2) 
r l 

T (4) 
r 4 

r (4) 

r 3 

(8) 


+ 
+ 

+ 

1 

+ 

+ 
+ 


V2 
V2 

v 3 

a 
a 2 
d 2 



Table 1: The operators and their representations. 



where the superscripts l,m mark different operators transforming according to the same 
representation. 

We will assume for the moment that there is only one representation of its kind. Let 
us then regard the full theory with quark and gluon operators first, before we go to the 
quenched approximation. In this case eq. (30) becomes 

( 0%) \_( Z™ Z™ \ ( Of (a) \ 

\ O^(fi) J \ Z'J'i Z'J'J J \ CP {a) J ' { ' 

where q (g) stands for quark (gluon). In perturbation theory one considers Green functions 
of the operators (5) - (8) computed between external off-shell quark and gluon states in 
a fixed gauge. We will impose the renormalization conditions 

(q(p)\0"^)\q(p)) = (q(p)\0"(a)\q(p)) tr 2 ee 2 , (32) 

(q(p)\0°(i,)\q(p)) = 0, (33) 
(g(p)\0«(fi)\g(p)) = 0, (34) 



{g(p)\0°(p)\g(p)) = (g(p)\Cy(a)\g(p)) 



tree 



(35) 



with \q(p)) (\g(p))) being a quark (gluon) state of momentum p. The understanding of 
eqs. (32) and (35) is that the conditions hold for those contributions which are propor- 
tional to the tree-level expressions. In the limit a — ► these conditions amount to the 
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continuum, momentum subtraction renormalization scheme. One obtains operators that 
are multiplicatively renormalizable by diagonalizing the anomalous part of the matrix of 
Z's in eq. (31). 

Let us now consider the quenched approximation. In this approximation 

(g(p)\O«(a)\g(p)) = 0. (36) 

This implies furthermore that Z q9 = 0. It then follows that for the quark operators we 
have 

O q (fi) = Z qq O q {a) (37) 

with the renormalization condition 

{q(p)\O q (ji)\q(p)} = (q(p)\O q (a)\q(p)) ^ = Z qq (q(p)\O q (a)\q(p)) . (38) 

Thus, in the absence of dynamical quark loops, the (singlet) quark operators do not mix 
with gluon operators under renormalization. 



4 Outline of the Calculation 

The task is now to compute (q(p)\0(a)\q(p)) for quark operators such as listed in table 1. 
General Considerations 

We use the standard Wilson action with r = 1. The calculation is done for zero renor- 
malized quark mass. We work in the Feynman gauge. The Feynman diagrams we have 
to calculate are listed in fig. 1. Only the operators which contain at least one covari- 
ant derivative receive contributions from the operator tadpole (fig. lb) and cockscomb 
(figs. lc,d) diagrams. 

The calculation basically amounts to the computation of integrals of the form 

d 4 k 



= J ■j—y [ IC lll ... lln (a,p,k), (39) 

where the integration is over the first Brillouin zone —it /a <k^< it /a. We follow Kawai 
et al. [7] and write 

1 = 1 +(1-1), (40) 

where 



d 



I(a,p) = I(a,0) + S Pa ^I(a, P ) p=o + ^g PaW ^-I(a, P ) 



1 ^ & 



2 



p=0 



1 d 3 

PoPpPiJ) — o — ^—Aa^P) +••• ( 41 ) 
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Figure 1: The lattice Feynman diagrams: (a) the vertex diagram, (b) the operator tadpole 
diagram, (c,d) the cockscomb diagrams, (e,f) the self-energy diagrams and (g,h) the leg 
tadpole diagrams. 
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is the Taylor expansion of the integral in the external momentum p. The order of the 
expansion is determined by the degree of ultraviolet divergence of X, so that the integral 
I — I is rendered ultraviolet finite. The latter can then be computed in the continuum, 
i.e. taking a — ► 0. The lattice integrals that remain to be computed do not depend on 
p, which greatly simplifies the calculation. However, the Taylor expansion about p = 
will in general create an infrared divergence. In order to regularize this divergence we 
use dimensional regularization with d > 4. The infrared singularities manifest themselves 
as poles in e = 4 — d. The infrared poles of X must cancel with those of I — I. In 
dimensional regularization X vanishes in the continuum limit. This can happen because 
the infrared divergence of X cancels exactly its ultraviolet divergence, thus turning the 
ultraviolet divergent integral X into an ultraviolet finite, but infrared divergent expression. 
The ultraviolet divergent contributions of the lattice integrals will cancel in the operators 
which we are interested in. 

The calculation is done analytically as far as this is possible. Using M athematica 
[12], we have developed a program package which, as input, only requires to state the 
Feynman diagram one wants to compute in symbolic form by specifying the operator, 
propagators and vertices defining the loop. The program then collects the various building 
blocks, works out the Dirac structure and does the necessary expansions. The continuum 
integrals are calculated analytically, including the infrared pole terms. The lattice integrals 
are separated into infrared singular and infrared finite integrals. The infrared singular 
integrals are computed analytically, while the finite integrals are computed numerically, 
as we will describe in the next section. The gamma matrices are chosen to anti-commute 
with 75, according to ref. [13]. This choice is not quite consistent with the axial anomaly 
[14]. But for the quenched calculation, as well as for non-singlet matrix elements, this 
should not matter. The results may differ though for the various schemes. Which scheme 
one will use finally will depend on the scheme that has been employed to compute the 
Wilson coefficients. 



Lattice Integrals 



The lattice integrals are brought by the program into the form 



r A d h t 2ni t 2ri2 t 2ri3 t 2rii 
B(n / ,n fl ;n 1 ,n 2 ,n 3 ,n 4 ) = a j { ^ )d gW/(o%w , (ojfc) , (42) 



where 





= 2 sin(afc |U /2) 


k 2 


= ^2 k % 






q(ak) 






A 1 


g(ak) 


= P. 



(43) 
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The integrals are infrared finite if 

n f + % < ni + ra 2 + ra 3 + ra 4 + 2. (44) 

If they are divergent, we proceed as follows to extract the infrared poles. We notice 
that 

^°^ + om - (45 > 

We then write 

^(n/,n a ;ni,n 2 ,n 3 ,n 4 ) = B(0, n f + n g ; n 1 , n 2 , n 3 , n 4 ) 

(46) 

+ [^(n/,n a ;ni,n 2 ,n 3 ,n 4 ) - #(0,n/ + ra 3 ; rai, ra 2 , ra 3 , ra 4 )]. 

The first integral, B(0,rif + ra 3 ; ni, ra 2 , ra 3 , ra 4 ), which is purely gluonic, can be computed 
analytically [15], giving 

K / n x 6(ra;rai,ra 2 ,ra 3 ,ra 4 ) ~ 

£5(0,ra;rai,ra 2 ,ra 3 ,ra 4 ) = h £5(0, n; n x , ra 2 , ra 3 , ra 4 ), (47) 

where 6 is some rational number divided by 7r 2 , and B can be reduced to a linear com- 
bination of one- dimensional elementary integrals, with exactly known coefficients, which 
can be computed to high precision. 

The degree of infrared divergence of the remaining integral, 

^ J J^y^ ^ ^ ^ (q n f(ak)g*9(ak) ~ g n f+ n °{ak) ' (4§) 

is reduced by two. By repeating this procedure several times, all infrared divergences are 
shoveled into purely gluonic integrals, which again can be computed analytically. This 
leaves us with infrared finite fermionic integrals of the form (42) only. 

As we already mentioned, these integrals will have to be computed numerically. Be- 
cause the final expressions involve quite a few of these integrals, in particular for the higher 
operators, we need to know their values to a high degree of accuracy. The first method we 
tried was integration by Monte Carlo using VEGAS [16]. This led to large errors, in par- 
ticular for the rank four operator. The second method we employed was four- dimensional 
Gauss-Legendre integration. Let us discuss and compare the two methods. 

For definiteness, let us consider the integral 

8(2, 2; 3, 0, 0, 0) = a* ^ (dfc g w ( 49) 

Because the integrand is an even function of k^, we may reduce the region of integration 
to the interval [0,7r/a]. Using 50 iterations with approximately 80,000 integration points, 
VEGAS gives the result 

B{2, 2; 3, 0, 0, 0) = 0.0037450(5). (50) 
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Figure 2: The lattice integral B(2, 2; 3, 0, 0, 0) as a function of 1/N 4 . The line is a linear 
fit to the data points. 

The Gauss-Legendre method distributes the integration points according to the zeroes 
of the Legendre polynomials. The density of integration points is particularly high near 
the endpoints. The endpoints themselves are, however, not covered. This feature is 
particularly well adapted to our type of integrals. We have computed the integral (49) 
using N = 30,35,40,45,50 and 55 integration points in each direction. We expect that 

B[N] = B[oo] + ^. (51) 

The result is shown in fig. 2. The agreement with our expectation is very good. We find 
similarly good results for all other integrals as well. We may therefore fit eq. (51) to the 
data and take B = B[oo] as our final result, with the error being given by its variance. We 
then obtain 

B(2, 2; 3, 0, 0, 0) = 0.0037455898(1). (52) 
This result is accurate to three more digits than the integration by Monte Carlo. 
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As an independent test we have computed the analytically known integral [15] 



5(0, 1; 1, 1,0, 0) = Z x = 0.1077813135399. 



(53) 



With our method we obtain 



Z t = 0.1077813135371(6), 



(54) 



which agrees with the result (53) up to the first 11 digits. 

5 Results 

In the quenched approximation and for the quark operators the renormalization constants 
can be written 



where 70 is the anomalous dimension of the operator. The same anomalous dimensions 
must appear in the corresponding Wilson coefficients, so that the product of Wilson co- 
efficient and operator matrix element is independent of 11. Obviously, the anomalous 
dimensions do not depend on the particular choice of representation of H(A) within a 
given 0(4) multiplet. 

Two Examples 

Before we state our results, let us present two examples which may serve to illustrate in 
what aspects the lattice calculation differs from the continuum calculation. 

The First Moment 

Let us first consider the operator 0^ v . To order g 2 we obtain 

i> 

{q(p)\0^\q(p)) = (1 + ci^Pj, + 027^ + c 3 7^ i/ p i/ + c 4 £8^ + c^p^ ~ traces, (56) 

where all c 4 - = 0(g 2 ). See table 8 in the appendix for values of the coefficients. Note 
that the contribution with factor C3 is non-0 (4) covariant. We shall look at two different 
representations: 



Z ((an) 2 ,g(a)) = 1 



16vr 2 



C F [7e>ln(a/i) + B a ] , C F = -, 



(55) 




This corresponds to the representation 73 '. We obtain from eq. (56) 
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1 h 

{q(p)\C{i4}\q(p)} = (1 + ci + c 2 )-(7ip 4 + 74Pi) + c 5 — Pip 4 - (57) 

z p 

The tree contribution, on the other hand, reads 

(<Kp)|C{i4}k(p))p e = ^( 7 iP4 + 74Pi)- (58) 

The standard renormalization procedure then leads to 

1 b 

Z (q(p)\£>{l4}\<l(p)} = Z (\ + Ci + C 2 )-(7lP4 + 74Pl) + Zc 5— PlPi 

z p 

1 h 
= ^(7iP4 + 74Pi) + c 5 — Pip 4 + 0(g 4 ). (59) 
z p 

Making use of the fact that Z = 1 + 0(g 2 ), we obtain 

Z = l- Cl -c 2 + 0(g 4 ). (60) 
There is no other operator transforming under Tg ' with which this operator could mix. 



044 " UOu + 22 + O; 



33 I 



This corresponds to the representation t[ 3 \ Here we have 

(q(p)\[0 44 - i(O n + 22 + 33 )]\q(p)) = (1 + ci + c 2 + c 3 )[74P4 - ^(TiPi + J2P2 



+73P3)] + c 5 [p^ - \(p\ +p 2 2 +P§)]4-(61) 



The tree contribution, on the other hand, reads 



1 



tree 



(q(p)\[0 44 - -(On + 22 + 33 )}\q(p)) = [ 74 p 4 - 3(71^1 + I2P2 + 73P 3 )]- (62) 
The standard renormalization procedure then gives us 

Z (q(p)\[0 44 - ^(di + C 22 + 33 )}\q(p)) = Z(l + ci + c 2 + c 3 )[74P4 - ^(liPi + I2P2 

+73P3)] + Zc 5 [pl - \{p\ + p\ + p^)]l 

6 p 

= [74P4 - ^(liPi + I2P2 + 73P3)] (63) 

+c 5 [pl- 1 -(pi+P 2 2 +pl)]4 + 0(g 4 ). 

p 

Making use of the fact that Z = 1 + 0(g 2 ), we obtain 

Z = 1 - ci - c 2 - c 3 + 0(# 4 ). (64) 
(3) 

There is no other operator transforming under ' with which this operator could mix. 

The origin of the difference of the renormalization constants (60), (64) is the presence 
of the non-0(4) covariant contribution to (56). 
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The Second Moment 



Let us now consider the operator Op, V p. To order g 2 we obtain 

{q{p)\0^p\q{p)) = (1 + cxYi^vVp + c 2 j 1/ Pf,p p + ctfpp^p,, + c^^vppl + c 5 jj w p 2 p 

+C6l P ^uPl + criS^nP^Pp + b^p-j^p^pv + b l ,p-{ v pp i p v ) (65) 

+C87fi7v7 P (pl + P 2 P ) + C9-^PfiPvPp H traces, 

where all c 4 - = 0(g 2 ). See table 9 in the appendix for values of the coefficients. The terms 
not listed here explicitly do not contribute to the specific operators which we will consider 
below. The contributions with factors C4 - c$ are non-0 (4) covariant. 

^{114} ~ l(^{224} + C{334}) 

(8) 

This corresponds to the representation . We obtain from eq. (65) 

(?(P)|[C{114} " ^(^{224} + C{334})]k(p)> = (1 + C X + C 2 + C 3 )^[27iPip 4 - J 2 P2P4 

-73P3P4 + l4Pl ~ 7}(l4Pl + 74Ps)] 

+(c 4 + C 5 + C 6 + 2c 8 )^[74p 4 

-^(74P2 + 74P3)] (66) 

+ C7[7lPlP4 - ^(72P2P4 + 73P3P4)] 



+c 9 ^p 4 [pl - ^(pI+pI)]. 
p z 



The tree contribution, on the other hand, reads 



1 



(?(P)|[C{114} - ^(C{224} + C{334})]k(p)> = ^llPlPi ~ I2P2PA ~ 73P3P4 



tree \ 



+74P? - ^(74P2 + 74P§)]- (67) 

This operator is not multiplicatively renormalizable. The non-0(4) covariant contributions 
with factors C4 — eg give rise to mixing with another operator of the same representation 
but of mixed symmetry. 

«411» ~ K C '«422)) + «433») 



The operator 0{ 114 y — ^((^{224} + ^{334}) mixes with the operator C^ 411 ^ — ^((^{{422)) 

C{{433))) which corresponds to the repre 
eq. (29).) From eqs. (29), (65) we obtain 

2( C '«422)) + <9«433»)]k(r-,/ l- ■ 2 



(8) 

C{{433))) which corresponds to the representation as well. {0^ vp ^ was defined in 
3qs. (29), (65) we obtain 

(?(p)|[<9«411» - ^(C{{422)) + ^((433)))]k(p)> = t 1 + c l ~ \{ c 2 + C 3 )](-27iPip 4 + 72P2P4 
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The tree contribution of this operator reads 

<?0>)|[<%11» " ^(<%22» + <%33»)]k(p)) 



+ 73P3P4 + 274Pi - 74P2 - 74P3) 

+ [c 4 - ^(c 5 + c 6 ) + 2c 8 ][2 7 4p? (68) 
-74P2 - 74Ps]- 



-271P1P4 + 72P2P4 + 73P3P4 

+2j 4 pj - j 4 pl - j 4 pl. (69) 



In an abbreviated form we express the renormalized operators as 
{} (ji) = Z {}{} O i} (a) + Z {}{{)) {{)) (a), 

= Z (0}{}°{}( a ) + Z «»«» c '«»( a )- ( 70 ) 

This is not a mixing in the usual sense as the matrix of anomalous dimensions is diagonal. 
The renormalization conditions that follow from (66) - (69) are 
11 3 1 

3 = Z {}{}^( 1 + C l + C 2 + C 3 + 2 C 7) " Z {}(())[ 1 + C l " + C3 )] 

11 3 
= Z {>{>3 + 3^ Cl + c 2 + c 3 + 2 C 7) - Z {}{{)) + °(9 4 ), (71) 

13 3 

= Z{}{}-(c 4 + c 5 + c 6 - -c 7 + 2c 8 ) + ^{}{{))[3 + 3ci - -(c 2 + c 3 ) 

+2c 4 - c 5 - c 6 + 4c 8 ] 
= ^(c 4 + c 5 + c 6 - ^c 7 + 2c 8 ) + 3Z {}(()) + 0(g 4 ), (72) 

where we have made use of the fact that = 0(g 2 ) and Z{}{} = l + 0(g 2 ). This then 

gives 



Z{}{} = 1 - ci - c 2 - c 3 - c 7 - -(c 4 + c 5 + c 6 + 2c 8 ) + 0(# 4 ), 
1 3 

Z {}{{)) = -q(c4 + c 5 + c 6 - -c 7 + 2c 8 ) + 0(g 4 ). (73) 



Similarly, one finds for and %((}}((}} 

^«»H = ~ 2C4 + C 5 + C 6 - 4C 8 + 0(# 4 ), 



%>«» = !- c i + 2 : ( c 2 + c 3 )-^(2c4-C5-c 6 + 4c 8 ) + 0(</ 4 ). (74) 



One can find other operators which do not give rise to mixing, but they require that 
the quark momentum is non-zero in more than one spatial direction. 
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The operator C^ 2 {i]4} 



The operator C^ 2 {i]4} requires special attention. To one loop it does not mix. However 
it turns out that (q(p)\0^ 2 ^ 4 y\q(p)} is linearly divergent, with the divergent part being 
given by 

2 

4.26568 -|— C F - 7 57 2 (74Pi - 7iP 4 )- (75) 
1671"^ a 

This contribution is directly proportional to r and vanishes for naive fermions. 

At first sight the occurrence of such a contribution might be surprising. But there is 
a good reason for that. In the OPE for §2 we have suppressed one operator, namely [17] 

°wZi]-»n} = "V (0 ^7[*7{ W ]75 B M • • • 5^} V - traces, (76) 

where m q is the quark mass. This explicitly quark mass dependent operator has also twist 
three. Usually, when the quark mass is renormalized multiplicatively, this operator can 
be neglected when one is only interested in the chiral limit. However, for Wilson fermions 
the situation is different. The divergent contribution (75) is exactly the contribution that 
is needed to renormalize the quark mass in the operator C^ 2 '{i]4}- 



Numerical Results 

We shall now present our numerical results. We will first consider the operators listed in 
table 1. In the appendix we shall also give results for (g(p)|(9 |UJ/ |g(p)), (q(p)\O l _ tl y P \q(p)) , 
(q(p)\0^ liyp \q(p)) and (q{p)\0 ^ p(J \q{p)) for general indices, from which the renormalization 
constants of all other representations [9] can be deduced. 

The results for the anomalous dimensions 70 are listed in table 2. The numbers agree 
with the anomalous dimensions known from the non-singlet Wilson coefficients [18]. 

The finite parts of the renormalization constants Bq are given in tables 3-6. Here 
we have also listed the individual contributions of the vertex, cockscomb, leg self-energy, 
leg tadpole and operator tadpole diagrams. (Remember that we are working in Feynman 
gauge.) The contribution of the leg tadpole diagram is 

8tt 2 £(0,1;0,0,0,0) = 8vr 2 Z 

= 12.23305015, (77) 

where we have used [15] Zq = 0.1549333902311. The total contribution of the operator 
tadpole diagrams is 

-n D 8ir 2 Z , (78) 

where ri£> is the number of covariant derivatives of the operator. There are ri£> operator 
tadpole diagrams. In the case of the operators 0^ 114 y — ^(^{224} + ^{334}) an( i ^{{411)) — 
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Operator 


lo 


C{14} 

C{44} - K {H} + ^{22} + C{33}) 
C{114} - |(C{224} + C{334}) 
«411» ~ K C '«422)) + «433») 
C{1144} + C{2233} - C{1133} ~ C{2244} 

o\ 

(I 5 

U [2{1]4} 


16 

3 

16 

3 

25 
3 

7 

3 

157 
15 



25 
3 

7 

3 



Table 2: The operators and their anomalous dimensions. 

l(C{{422)) + C{{433))) the operator tadpole contribution is spread among and Bs^^y 

Here the result (78) holds for ci, the order g 2 contribution with tree structure, and hence 
for — and -£?{{)){{)) — {{)){}• The °P era t° r tadpole diagrams have the 

opposite sign to the leg tadpole diagrams. In the case of the operator O^, which contains 
one covariant derivative, leg tadpole and operator tadpole diagrams cancel exactly. In all 
other cases the tadpole diagrams account for more than 60% of the total contribution. 

The renormalization constants of the operators C{i4} and — ^(O^ny + Os^2} + 

^{33}) differ only in the vertex contribution to the finite part. 



Conversion to the MS Scheme 



The Wilson coefficients are usually computed in the MS scheme, so that one would need 
to know the renormalization constants in this scheme too. The result in the MS scheme is 
easily obtained. In table 7 we give the finite contribution of the continuum integrals to Bq, 
which we denote by Bq 71 . Here 7^ = 0.57721566 is Euler's constant. The renormalization 
constants in the MS scheme are then given by 

B^ s = B Q -B c S n . (79) 

Similarly, the renormalization constants in the MS scheme are obtained by 

Bp = B - B c o 0n + ^( 1E - In Air). (80) 
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Diagram 




^C{44}-|( C '{ll}+ C '{22}+ C '{33}) 


Vertex 
Cockscomb 
Leg Self-Energy 

Leg Tadpole 
Operator Tadpole 


2.2930524(2) 
-5.0772671(1) 
-0.3806456(7) 
8tt 2 Z 
-8tt 2 Z 


3.5753197(3) 
-5.0772671(1) 
-0.3806456(7) 
8tt 2 Z 
-8tt 2 Z 


Total 


-3.1648603(2) 


-1.8825929(3) 



Table 3: The operators 0{i4} and (9{ 44 } — |((9_r n } + 0{ 22 } + ^{33})- 



Diagram 


B {}{\ 


%«» 


%>{} 


%>«» 


Vertex 
Cockscomb 
Leg Self-Energy 
Leg Tadpole 
Operator Tadpole 


1.357071(1) 
-6.73958572(2) 
-0.3806456(7) 
8tt 2 Z 

2^2 64^2^ 


-0.0559027(9) 
0.21263441(7) 

2^2 16^2^ 


-1.2849696(4) 
7.4327048(3) 

4vr 2 - 32vr 2 Z 


1.0080635(2) 
-2.0391253(2) 
-0.3806456(7) 
8vr 2 Z 

4^2 80 ^2 y 


Total 


-19.571840(1) 


-0.36847847(9) 


-3.3060478(2) 


-16.7960186(1) 



Table 4: The operators {114} - |(C{ 2 24} + C{334}) and 0« 4 n)) - |(C{{422)) + «433»)- 



Diagram 


-^^{1144} +^{2233} -^{1133} -^{2244} 


Vertex 
Cockscomb 
Leg Self-Energy 

Leg Tadpole 
Operator Tadpole 


-6.37131(2) 
-5.860427(1) 
-0.3806456(7) 
8tt 2 Zo 
-24tt 2 Zo 


Total 


-37.07849(2) 



Table 5: The operator C {1144} + C {2233} - C {1133} - C {2244} . 
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Diagram 

o 




^{214} 


r2-rn4> 


Vertex 
Cockscomb 
Leg Self-Energy 

Leg Tadpole 
Operator Tadpole 


3.94387868(6) 

-0.3806456(7) 
8tt 2 Z 


0.661669(1) 
-7.6095645(2) 
-0.3806456(7) 
8tt 2 Z 
-16tt 2 Z 


0.9887912(2) 
-4.0525425(2) 
-0.3806456(7) 
8tt 2 Z 
-16tt 2 Z 


Total 


15.79628324(6) 


-19.561590(1) 


-15.6774470(1) 



Table 6: The operators Q\, C^ 2 i4} an< ^ ^2{i]4} - 



Operator 



TDCOTl 

a 



o 



{14} 



C{44} - K {H} + ^{22} + C>{33}) 



{114} 



o 



«4H» 



■ K C '{224} + C{334}) 
|( C '«422)) + «433») 



{1144} 



o 



{2233} 



^{1133} ~~ ^{2244} 



(I 5 

^{214} 
U [2{1]4} 



40 
" 9 
40 
" 9 
67 
" 9 

35 
" 18 
2216 
" 225 



f + §( 7 £-m47r) 
§( 7s - ln47r) 
f ( 7s -ln47r) 
f| + f( 7 £-m47r) 
W(7s-ln4vr) 



9 

35 
" 18 







-f + f( 7s -ln4vr) 
|( 7S - ln47r) 



Table 7: The finite contribution of the continuum integrals to Bq. 
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6 Summary 



We have computed the renormalization constants of the leading twist lattice bilinear quark 
operators up to spin four. The calculation was done in the quenched approximation using 
Wilson fermions with r = 1. Results for other values of r can be obtained from the authors. 

For non-singlet quark operators and to one-loop order there is no difference between the 
quenched approximation and the full theory including dynamical quarks. The difference 
shows in the singlet operators, and here only in the renormalization constants Z qg and 
Z'J'J in eq. (31). 

The renormalization constants for the axial vector current and for the operator 
^{14} were known before [19, 20]. The calculation of the renormalization constant for the 
operator C{n4} — |(C{224} + ^{334}) and its mixing parameters was done in parallel [21] 
to ours [2]. These authors use a slightly different basis of operators though. The results 
all agree. 

We have explicitly stated the contributions that come from the tadpole diagrams. It 
is then straightforward to compute the renormalization constants in tadpole improved 
perturbation theory [3]. One possibility is to write 

1 " ^ CfB ° = un F~ X {l ~ 4~2 CfB ^ = " TT-2 C ^o) + 0(g* 4 ) (81) 

(np>: number of covariant derivatives), where 

uf = ^TrE^ = 1 - j^C F 8tt 2 Z + 0(g 4 ), (82) 
Ufj, being the link matrix in Feynman gauge, and 

B = B + (n D - 1)8tt 2 Z . (83) 

This reflects the fact that one finds ri£> operator tadpole and one leg tadpole diagrams, 
which are of the same magnitude but have opposite sign. In the case of the operators 
C{ii4} - U°{224} + C{334}) and C«4ii)) - ^(C«422)) + C«433))), which mix, one has to 
consider — and -£?{{)){{)) — instead. One factor of up will be absorbed 

into the normalization of the quark states, k c — ► k c = K c up, because k c is expected to have 
a better behaved perturbation series, k c = |(1 + ^-0.066) + 0(g* 4 ). As the expansion 
parameter g* one uses the coupling constant renormalized at some physical scale. In ref. [5] 
we have compared our results with tadpole improved perturbation theory. 
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Appendix 



In this appendix we shall present our results for the full tensors. We are only interested 
in the finite contributions. We write 

g 2 

(q(p)\0^(a)\q(p)} = 7 ^ + -—^C F -traces, (84) 

g 2 

(q(p)\0^p(a)\q(p)} = -)^p v p p + -—zCfB^p - traces, (85) 



g 2 

{q{p)\Ol vp {a)\q{p)) p2=1/a2 = l»l5PvP P + J^2 C F B lv P ~ traces, (86) 



g 2 

{q{p)\0^ pa {a)\q{p)) = lp,p v p p Pa + -r^C F B^ pa - traces, (87) 

where the logarithmic (anomalous) contributions are switched off, and 

B ^ = (88) 

Blp = ECp c5 ' (9°) 

Bjifpa — ^ ^ i jj,fpa (^1) 

The numerical values of the coefficients c, c 5 for the finite contributions are listed in tables 
8 - 11. We have omitted the singular contributions of order 1/a 3 (in B^ vp(T )^ 1/a 2 (in 
up an( i B^vp) an d order 1/a (in Bp Vpa , B^p, B^ vp and Bp V ). They will cancel 
out after symmetrization and subtraction of the traces. The singular contributions may 
be obtained from the authors. 



tpll> 


c 


IP^V 


2.0523298(1) 




1.2822674(5) 


lv Pn 


-1.5824302(1) 


111 Pv 


-1.5824302(1) 


/PPnPv/p 2 


-4/3 



Table 8: The operator O l 
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C 




c 


Kp ip Pp 2 


0.493942(1) 




-0.00564(1) 


Kp lu Pp 2 


0.138023(1) 


K» IpPp 2 


-0.886064(1) 


IpSppPv 


1.8345574(6) 


Pip IpPu 


-0.43037257(4) 


Kp IpPp Pu 


0.360814(2) 


Kp lu Pp Pu 


0.360814(2) 


IpPp Pu 


-0.6741686(7) 


K P ip Pu 2 


-3.051131(1) 


Kp lu Pu 2 


0.493942(1) 


lp lu IpP 2 


0.5120432(1) 


IpSpvPp 


0.9738123(5) 


Pip lu Pp 


0.43037257(4) 


K" lp Pp Pp 


0.360814(2) 


luPp Pp 


-1.5349138(6) 


IpPu Pp 


-16.7985438(6) 


Kv IpPp 2 


-0.5301448(9) 


2 

lp lu Ip Pp 


n k 1 on a oo/ 1 \ 
U. 0lzU4oz^l ) 


PPp Pu Pp/P 

ip <W p 2 


-1 


Kp ipP 2 


-0.3175361(2) 


-0.0614063(8) 


Kp IvP 2 


0.4422501(3) 


K» IpP 2 


-0.2590425(2) 


lp lu Ip P 2 


0.35064632(4) 


/pKpPp 


0.2185109(6) 


/PPp bp V p 


-0.494257(2) 







Table 9: The operator O \ 



t 5 

l pup 


c 5 


t 5 

l pup 


c 5 


i>K P is Pp 


-0.1618359(5) 


£>15 Pp K»P 
lp 15 Pp 2 <W 


-0.329227(2) 


Kp lp 15 Pp 2 


0.064073(1) 


0.14207(1) 


Kp Iv 75 Pp 2 


0.29828(1) 


K» 7p 75 Pp 2 


-0.725806(1) 


t>Kp 75 Pu 


1.644384(5) 


Pip lp 15 Pu 


-0.43037257(4) 


Kp Ip. 75 Pp Pu 


0.040298(2) 


Kp lu 75 Pp Pu 


0.525844(2) 


lp 15 Pp Pu 


-0.864342(6) 


Kp lp 75 Pu 2 


-3.428208(1) 


Kp lvl5 Pu 2 


0.654199(1) 


lp lu lp 15 Pu 2 


0.5120432(1) 


i>Kv 75 Pp 


0.7836389(5) 


Pip lu 15 Pp 


0.43037257(4) 


K" lp 75 Pp Pp 


0.040298(2) 


lu 15 Pp Pp 


-1.7250872(6) 


lp 15 Pu Pp 


-16.9721616(5) 


K» Ip 75 Pp 2 


-0.369887(9) 


lp lu lp 15 Pp 2 


0.5120432(1) 


Pl5 PpPu Pp/P 2 
lp 15 <W P 2 


-1 


Kp ip i5 p 2 


-0.1218283(2) 


-0.160058(1) 


Kp lv 75 P 2 


0.6324235(3) 


K" 7p 75 P 2 


-0.0688692(2) 


lp lu lp 15 P 2 


0.35064632(4) 







Table 10: The operator O b 
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C 




C 


phpPv Pv 


-3.08982 


1) 


t>Pa Pv hpp 


15.10485 


[4) 


hp 7(7 lp P 


-0.2228411 


v 4 ) 


t> hp la lp Pa Pv 


-1.4585045 


;5) 


P hp lp lv Pa Pv 


-1.4585045 


;5) 


c 9 
hp 7(7 Pa Pv 


6.77201( 


2) 


7(7 Pa Pv hpp 


-76.2628 


i 2 ) 


c 9 

hplp Pa Pv 


7.5907( 


2) 


c c 2 
hphplvPa Pv 


-26.80781 


;6) 


c 9 
hplp Pa Pv 


7.35454( 


2) 


c 9 

<Vp 7(7 lp lv Pa Pv 


4.974603 


i 2 ) 


hp lp lv P 


-0.2228411 


[4) 


O 

7(7 lp IpPa Pv 


-0.4352333 


v 4 ) 


phpPpPv 


-3.84512( 


1) 


hp 7(7 Pa Pp Pv 


11.55906 


[A) 


hp lp-P<y Pp. Pv 


15.55123 


[4) 


hp lv Pa Pp Pv 


14.47606 


[4) 


IpPa Pp Pv 


-4.66978( 


1) 


c 9 
hp 7(7 Pp Pv 


5.45522 


[2) 


c 9 

hplpPp Pv 


7.90623( 


2) 


J c c 2 
Php hp Pv 


7.73053 


[2) 


/ c 9 
Php la IpPv 


-0.699931 


;i) 


/ c 9 
Phplp lv Pv 


-0.5068522 


V) 


c c 2 
hp hp la Pa Pv 


-23.70894 


;6) 


c 9 
hp lp Pa Pv 


7.46065 


[2) 


c 9 

hp lv Pa Pv 


5.90631( 


2) 


9 c 

7^ P<7 Pv hpp 


-26.93646 


;6) 


c 9 
hplalp lv Pa Pv 


4.866452 


[2) 


o 

7a* 7^ IpPa Pv 


-0.4352333 


V) 


c 9 

hp la Pp Pv 


4.9362( 


2) 


c 9 
Cap lv Pp Pv 


6.92001 


[2) 


hp 7(7 Pv 


1.730699 


;6) 




-7.89989 


[2) 


hplp Pv 


1.134569 


[7) 


hp hp lv Pv 3 


-28.7742 


;6) 


hplp Pv 


2.557857 


;6) 


hp 7(7 7a* lv Pv 


3.764515 


;5) 


7(7 lp lp Pv" 


-0.0321032 


[4) 


i>hv Pa Pp 


-4.39407 


;i) 


t>Pa Pp hpv 


15.11438 


[4) 


t>lp lv Pa Pp 

9 c 

7(7 Pa Pp <Wj, 

c 9 

<W 7^ Pa Pp 


1.0825005 


V) 


c 9 

<V7<7 Pa Pp 


8.31511( 


2) 


-76.258 


[2) 


c 9 
hvlp Pa Pp 


6.3978( 


2) 


7.35454 


[2) 


9 

7(7 lp lv Pa Pp 


-1.317318 


;i) 


fchvPpPp 


-3.84512 


;i) 


hv la Pa Pp Pp 


11.54952 


[4) 


hv IpPa Pp Pp 


15.55123 


[4) 


lv Pa Pp Pp 


-5.53052( 


1) 


c 9 

hv la Pp Pp 


7.35454 


[2) 


c 9 

hvlp Pp Pp 


7.90623( 


2) 


phpPv Pp 


-4.39407 


;i) 


Pla IpPv Pp 


1.0825005 


;6) 


hp la Pa Pv Pp 


15.55123 


[4) 


IpPa Pv Pp 


-25.11451 


;i) 


la Pp Pv Pp 


-3.36552 


;i) 


c 9 
hplv Pv Pp 


8.31511( 


2) 


9 

7(7 lp lv Pv Pp 


-1.206783 


;i) 


c 9 
hv IpPa Pp 


8.31987( 


2) 


lp Pa Pp 2 hpv 


-26.95077 


;6) 


lp lv IpPa Pp 2 


-u.zoyou i 


i j 


h,, In V,, Z>„ 2 

^ av I p r p r p 


6.92001 


[2) 


hn In Vi, Z>„ 2 

a p ] p rv r p 


8.31987( 


2) 


la lp IpPv Pp 2 


-0.259507 


;i) 


hv la Pp 3 


2.557857 


;6) 




-8.12096 


[2) 


hvlp Pp 3 


2.707691 


[7) 


hp lv Pp 3 


2.557857 


\<o) 


la lp lv Pp 3 


-0.4456823 


;9) 


t>Pa PpPv Pplp 2 


-4/5 


t>hp hv P 2 


-2.656007 


;5) 



Continued on next page 
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/ 
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D.OOUD ( 
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/ 
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( .0440 ( 


9^1 

/ 


71 — S sr i ii i r\ 
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70 n7Q8 

- , o.u , yo 


Z J 
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-U.0U44DOO 


o J 


h 7t 7„ VrP 8n,,n 

/¥ 1 (7 1 pi ¥<7 ^ (7 is p 


-u.uuzoooyo 


1 1 


h S„„ 7„ V J 2 

It ^ a p lp if L J a 


-u.oyyyoi 


1 ^1 

1 1 


8rrn 8 „u In Vn^ 

(7 p p,U j (7 ¥ <7 


98 77zL9 
-Zo. , , ^rZ 


J 


fan Snn la Va^ 

(7 U p,p /(J 1 (J 


91 9zL^1 7^ 
-Z 1 .Z40 1 ( 


/ 


8nn 8,,n In Vn ^ 
^ (7 pi ^ is p j (7 ¥(7 


-Z0.4DU1D 


D J 


'y ii ^ 8 

1(7 ¥(7 (7 pils p 


CQ/i 777/ 
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1 \ 
1 1 


8 ^ 8 , , i , ^) 

^(7 p ^ pils /¥ 


9 KK87^/I 
Z.DDo , 04 


D J 


8,,n 7m Vn^ 

^ is p 1 pi ¥(7 


z.ooooy4 


8^1 
o ) 


P> P<7 8n\jp 


"1U.OUZZU 


'8"l 


""7 pi P(7 8 p,is p 


7 8QQ8Q/ 
- , .oyyoy 


) 


8,,n lis Vn 3 

p,p j U ¥ <7 


1 KKK/IQ9 

i .ooo^yz 


D J 


7 71 ^ 8 

1 f ¥ (7 (7 pip 


9/1 /iQ/in^/ 

-Z4.4y4U0 




8 hi, 7 n Vn^ 
^ piu j p ¥<7 


9 ^^78^7 


D J 


""7 ' p P(7 8 (j 


9zL 97908/ 
-Z4.z j zyo 


7\ 


8nn In In 7iy Vn ^ 
^ (7 p 1 (7 1 pi 1 Is ¥<7 


^ 87^flzLQ 

o.o , ou^y 
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In 7iy In Vn 
1 pi 1 is 1 p ¥(7 


-u.uoziuoz 


A\ 
4 1 


8 n 1/ 8 n n 

(7 U p,p /¥ 


2.223071 


;6) 


if 8,,n Vn Vn 
/¥ u p ¥ (7 ¥ pi 


-3.84512( 


1) 


IPPa Pp hvp 


12.65199 


[A) 


t> Pa Pp Spvp 


15.10485( 


4) 


Sap Sup jp 


2.668754 


\<o) 


hp la Pa 2 Pp 


5.90631( 


2) 


la Pa 2 Pp, t>avp 


-63.6107 


[2) 


la Pa 2 Pp Spup 


-27.148( 


6) 


flap fipvlp Pa 2 Pp 


-23.71333 


\<o) 


SauSpp lp Pa 2 Pp 


-23.71809( 


6) 


flap lu Pa 2 Pp 


6.33606 


[2) 


hulp Pa 2 Pp 


6.33606( 


2) 


fpp la lp lv Pa 2 Pp 


0.2210686 


[2) 


Sapvp IP 


-8.90574( 


2) 
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tpvpa 


C 


tpvpa 


C 


fabrrn S ,,,, V , , ^ 

lr a p pv i j p 


7.34437(2) 


h 8m, S ,, „ v,?' 

Ir av pp I p 


6.33067(2) 


^ (7 p ^ pv ] (7 l J (T l J p 


-26.80304(6) 


Sm, S,,n 'Yrr V rr V ,, ^ 

^av pp la i J a i J p 


-17.05583(6) 


Kp lp Pa Pp 2 


7.90623(2) 


IpPa Pp 2 fiavp 


-23.85046(6) 


lp Pa Pp? bpvp 


-76.2672(2) 


Spp lv Pa Pp 2 


4.9362(2) 


bpv IpPa Pp 2 


7.34977(2) 


hp la Pp 3 


1.666492(6) 




-5.68185(2) 




-24.27298(7) 


flap fipv lp Pp 3 


-25.46016(6) 


Sav &pp lp Pp 3 


-25.68122(6) 


Saplv Pp 3 


2.112174(6) 


t>avlp Pp 3 


2.112174(6) 



Table 11: The operator 0^ vpcT . Here jf) 3 = J^^lpP^ 
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